均匀磁场下碳离子笔形束的剂量变化分析及位置修正方法

分析外加均匀磁场对于碳离子笔形束剂量分布的影响,并考虑修正这种影响,为磁共振成像引导碳离子放射治疗的临床应用提供指导。本文利用蒙特卡罗方法模拟计算了不同能量碳离子笔形束在不同强度磁场下的剂量分布情况,发现垂直于碳离子束入射方向的均匀磁场对于碳离子笔形束射程缩短的影响很小,磁场对碳离子束的主要影响是引起束流横向偏转,特别是碳离子束布拉格峰位置的横向侧移。横向侧移程度与碳离子束的能量和磁场强度相关,根据模拟结果,得到了一个计算碳离子束布拉格峰在磁场中相对横向偏转的方程,并提出一种校正外加磁场引起的碳离子束布拉格峰横移的角度修正方法。这些结果可用于指导磁共振图像引导碳离子放射治疗计划系统的研发。     

80.5 MeV/u碳离子诱发铜靶的放射性剩余产物测量

高能重带电粒子能直接穿透靶原子核外电子层,与原子核发生直接碰撞,发生散裂反应,产生一系列具有放射性的剩余产物核.重带电粒子诱发靶材放射性剩余核与辐射防护和人员安全有着密切联系,当前,大部分剩余核产额主要依靠蒙特卡罗粒子输运程序进行模拟计算,其准确程度亟需通过实验测量进行准确评估.本文利用能量为80.5 MeV/u的(12)^C⁶⁺粒子对薄铜靶开展了辐照实验与伽玛射线测量,结合伽玛谱学分析方法,得出了辐照产生的18种放射性剩余产物的初始活度和产生截面值,并与PHITS模拟结果进行对比.结果表明,PHITS模拟程序对放射性剩余核种类的估计具有较高可靠性,在其绝对产额方面,与实验测量仍具有较大偏差.     

Monte Carlo simulation of electric conductivity for pure and doping fullerene (C60)

In this paper, we have used Monte Carlo (MC) method to simulate and study the temperature and doping effects on the electric conductivity of fullerene (C60). The results show that the band gap has reduced by the doping and the charge carrier transport is facilitated from valence band to conduction band by the temperature where is touched a 300 K. In this case, the conductivity reached a value of $4\times {10}^{-7}\text{S}{\text{cm}}^{-1}$. The electric conductivity of C60 can increase by the triphenylmethane dye crystal violet (CV) alkali metal to reach $4\times {10}^{-3}\text{S}{\text{cm}}^{-1}$ at 303 K. Our results of MC simulation have a good agreement with those extracted from literature [10], [33].     

Monte Carlo study of interfacial silicon suboxide layers and oxidation kinetics

A simple simulation scheme that simultaneously describes the growth kinetics of SiO₂ films at the nanometer scale and the SiO_x/Si interface dynamics (its extent, and spatial/temporal evolution) is presented. The simulation successfully applies to experimental data in the region above and below 10 nm, reproduces the Deal and Grove linear-parabolic law and the oxide growth rate enhancement in the very thin film regime (the so-called anomalous region). According to the simulation, the oxidation is governed mainly by two processes: (a) the formation of a transition suboxide layer and (b) its subsequent drift towards the silicon bulk. We found that it is the superposition of these two processes that produces the crossover from the anomalous oxidation region behavior to the linear-parabolic law.     

锥形光纤中光波传输特性的蒙特卡罗模拟

为了分析锥形光纤中光子流的传输特性和空间分布特性，采用了蒙特卡罗算法对其进行了模拟仿真，得到了大量光子经过锥形光纤后的光子分布图.模拟实验结果表明，通过锥形光纤的大量光子集中在探测器中心，到达探测器某一位置的光子数随着该位置到探测器中心距离的增加而减少，并且在锥形光纤锥度比一定的条件下，改变锥形光纤长度，光强度随着长度的增加而下降，但通过锥形光纤的光子分布是一样的.这些研究结果对锥形光纤的设计和实际应用提供了重要依据.     

Non-universal critical behaviour of heteronuclear dimers on a square lattice––A Monte Carlo study

Monte Carlo simulation within the grand canonical ensemble, the histogram reweighting technique, and finite size scaling analysis are used to explore the phase behaviour of heteronuclear dimers, composed of A and B type atoms, on a square lattice. We have found that for the models with attractive BB and AB nearest-neighbour energy, u_BB=u_AB=−1, and for non-repulsive energy between AA nearest-neighbour sites, u_AA<0, the system belongs to the universality class of the two-dimensional Ising model. However, when u_AA>0, the system exhibits a non-universal critical behaviour. We have evaluated the dependences of the critical point characteristics on the value of u_AA.     

电子回旋共振放电的电离特性PIC/MCC模拟(Ⅰ)——物理模型与理论方法

采用粒子模拟与蒙特卡罗相结合(PIC/MCC)的方法对电子回旋共振(ECR)放电中的电离过程进行了模拟，其中带电粒子与微波的相互作用由PIC方法的电磁模型描述，粒子间的碰撞过程由MCC方法描述.考虑的碰撞类型有电子与中性粒子的弹性、激发、电离碰撞，离子与中性粒子的弹性、电荷交换碰撞，碰撞截面均依赖于能量而变化.阐述了理论分析的过程，为数值模拟ECR放电奠定了基础. 关键词： 电子回旋共振放电 粒子模拟 蒙特卡罗 电离     

Approximation of photonuclear interaction cross-sections

Photonuclear interaction cross-sections from the GEANT4 database are approximated for all nuclei and all energies (from the hadron production threshold to about 40 TeV). The approximation methods in the giant-dipole resonance region, nucleon resonance region, and high-energy region are improved with respect to existing approximations. As an application of the approximation for photonuclear cross-sections, an improved method of calculating electronuclear cross-sections is developed. The interaction cross-section of virtual photons with nuclei at high Q² are approximated and a simple algorithm for describing the electronuclear reactions, including high-Q² scattering, is proposed. Received: 22 February 2002 / Accepted: 6 May 2002     

一种高信噪比的双光束激光光声光谱仪原理及应用

作为一种新兴的检测手段 ,激光光声光谱技术与其他检测技术相比具有很多优点。本文设计的光声光谱仪用激光做光源 ,有两个光声池 ,分别用来放置参考样品和待测样品 ,输出结果为二者光声信号的比值。利用这一比值 ,可以由参考样品的性能参数方便地求出待测样品的相关性能参数。该光声光谱仪有效地减小了本底吸收噪声的影响 ,提高了信噪比 ,扩大了固体光声理论的应用范围。本文还阐述了该光声光谱仪在定性定量分析中的一些应用     

Construction algorithms for polynomial lattice rules for multivariate integration

We introduce a new construction algorithm for digital nets for integration in certain weighted tensor product Hilbert spaces. The first weighted Hilbert space we consider is based on Walsh functions. Dick and Pillichshammer calculated the worst-case error for integration using digital nets for this space. Here we extend this result to a special construction method for digital nets based on polynomials over finite fields. This result allows us to find polynomials which yield a small worst-case error by computer search. We prove an upper bound on the worst-case error for digital nets obtained by such a search algorithm which shows that the convergence rate is best possible and that strong tractability holds under some condition on the weights.

We extend the results for the weighted Hilbert space based on Walsh functions to weighted Sobolev spaces. In this case we use randomly digitally shifted digital nets. The construction principle is the same as before, only the worst-case error is slightly different. Again digital nets obtained from our search algorithm yield a worst-case error achieving the optimal rate of convergence and as before strong tractability holds under some condition on the weights. These results show that such a construction of digital nets yields the until now best known results of this kind and that our construction methods are comparable to the construction methods known for lattice rules.

We conclude the article with numerical results comparing the expected worst-case error for randomly digitally shifted digital nets with those for randomly shifted lattice rules.